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Abstract 

A unifying scheme based on an ancestor model is proposed for generating a wide 
range of integrable discrete and continuum as well as inhomogeneous and hybrid models. 
They include in particular discrete versions of sine-Gordon, Landau-Lifshitz, nonlinear 
Schrodinger (NLS) , derivative NLS equations, Liouville model, (non-)relativistic Toda 
chain, Ablowitz-Ladik model etc. Our scheme introduces the possibility of building a novel 
class of integrable hybrid systems including multi-component models like massive Thirring, 
discrete self trapping, two-mode derivative NLS by combining different descendant models. 
We also construct inhomogeneous systems like Gaudin model including new ones like 
variable mass sine-Gordon, variable coefficient NLS, Ablowitz-Ladik, Toda chains etc. 
keeping their flows isospectral, as opposed to the standard approach. All our models are 
generated from the same ancestor Lax operator (or its q 1 limit) and satisfy the classical 
Yang-Baxter equation sharing the same r-matrix. This reveals an inherent universality in 
these diverse systems, which become explicit at their action-angle level. 
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1. INTRODUCTION 



Though integrable models represent only a special class of nonlinear systems, their numbers 
and varieties discovered till today have become amazingly large. Therefore it is particularly 
important now to have well defined schemes, which will be able to generate them in a system- 
atic way, find out their interrelations, detect the fundamental ones and identify their universal 
properties. Reduction of Lax operators in AKNS spectral problem [1], classification of soli- 
ton bearing equations through self-dual Yang-Mills equation [2] , gauge unification of nonlinear 
Schrodinger (NLS)-type models [3] are few of such successful approaches. However most of 
these schemes are designed to deal with the continuous models only, whereas the importance 
and significance of discrete integrable systems have been well emphasized in recent years [4]. 
Moreover, the algebraic approach in classical integrable models, though has a rich and sophis- 
ticated formulation through the classical Yang-Baxter equation and the classical r-matrix [5], 
as it appears, has not been exploited fully. 

Our aim here is therefore to propose an unified algebraic scheme for systematic generation 
of a large class of integrable discrete models, based on their underlying Poisson bracket (PB) 
structure. The specialty of this class of models is that , they can be easily quantized to yield 
the corresponding quantum integrable systems and their classification may be done through the 
associated classical r-matrix with its known trigonometric and rational solutions. We present 
an integrable discrete ancestor model linked with the trigonometric r-matrix (and its — > 1 
form related naturally to the rational solution of r) and containing a set of arbitrary parame- 
ters. Various choices of these external parameters define in turn different underlying algebraic 
structures and the associated Lax operators. This generates through suitable realizations a 
wide range of diverse integrable systems sharing the same r-matrix with their ancestor model. 
They are by construction integrable discrete models with few of them having also well defined 
field limits. 

Our scheme, the basic idea of which is borrowed from the quantum domain [6, 7], appears 
to be effective not only in classifying an important class of discrete models as well as their field 
limits, but also their inhomogeneous extensions. Along with the exactly integrable discrete 
versions of the well known models like sine-Gordon, Landau-Lifshitz equation, NLS, deriva- 
tive NLS (DNLS), Liouville model, relativistic and nonrelativistic Toda chain, Ablowitz-Ladik 
model, we also obtain new inhomogeneous models like variable mass sine-Gordon and more gen- 
eral variable coefficient NLSs, as well as Gaudin model, inhomogeneous Ablowitz-Ladik model 
and Toda chains. As an important application of our scheme we may construct novel families 
of integrable hybrid models, by combining different descendant models in different domains of 
the lattice space ^, or by fusing copies of a single component model to get its multi-component 
generalization. Moreover, the present method of generating integrable inhomogeneous discrete 

^Very recently such field models attracted attention [8] 
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and continuum models reveals the intriguing fact that the conventional approach by considering 
space-time dependent spectral parameter X{x,t) [9]- [14] is rather restricted and even appears 
to be misleading, since it would lead in general to a dynamical r-matrix spoiling the underlying 
algebraic structure and forbidding therefore the possible quantization of the models and their 
usual action-angle formulation. Moreover, for more general inhomogeneous sine-Gordon and 
NLS models, as we find here, the conventional treatment of nonisospectral flow would likely to 
fail. In our approach on the other hand the necessary isospectrality is kept intact by taking 
constant A as in the original homogeneous case and the inhomogeneity is introduced through 
arbitrary parameters, which act like Casimir operators in the associated Poisson algebra. 

Since all these models, in spite of their manifestly diverse forms and nature, are generated 
from the same ancestor model sharing the same r-matrix (or its g ^ 1 form), it reveals an in- 
triguing universality among them which is reflected prominently in their description of complete 
integrability through action-angle variables. 

The paper is arranged as follows. In sec. II we review the theory of integrable systems 
satisfying classical Yang-Baxter equation associated with classical r(A — fi) matrix. Sec. Ill 
presents the explicit form of the ancestor model and its g ^ 1 limit together with the underlying 
PB algebras. We introduce our generating scheme in sec. IV and construct concrete models. 
Sec. V accounts for the generation of integrable inhomogeneous as well as hybrid models. Sec. 
VI focuses on the universal property of all descendant models by explicit construction of their 
action-angle variables. Sec. VII is the concluding section. 

11. CLASSICAL YANG-BAXTER EQUATION AND INTEGRABLE SYSTEMS 

By integrabihty of a nonhnear discrete system defined on a lattice with sites j — 1,2,. . .,N, 
we mean it in the Liouville sense by requiring the existence of its number of independent 
conserved quantities Cn, n — 1,2, ... ,N including the Hamiltonian H of the system with the 
criteria {Cn, Cm} — 0. Such conserved quantities can be considered as the action variables 
generated from a spectral parameter A-dependent transfer matrix as: r(A) = J2n=i C'nA'* and 
consequently the integrability criteria may be replaced by the single condition 



For deriving this condition therefore along with the conventional linear spectral problem: 
Tfe+i(A) = Lfc(A)Tfc(A) we define also the PB algebra for its Lax operator Lk{X) in a specific 
form, which is known as the classical Yang-Baxter equation (GYBE) [5] 



associated with the classical r(A — /i)-matrix playing the role of structure constants. For the 
associativity of algebra (2) ensuring its Jacobi identity, the r-matrix in turn must satisfy another 



{r(A),r(/x)} = 0. 



(1) 



{Lfc(A)(g), Liifi)} = Ski[r{X - fx), Lfc(A) O Lfc(/x)] 



(2) 
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form of CYBE: 



[ri2(A - //), ri3(A - S)] + [r^iX - fJ,),r23{f^ - S)] + [ri3(A - S),r23{f^ - S)] = 0. 



(3) 



It is crucial to observe that, though there is a variety of Lax operator solutions to (2) with 
different basic operators and spectral parameter dependence, representing a wide range of 
integrable systems (for a list see sec. IV), the associated r-matrix solutions satisfying (3) are 
only of three types: elliptic, trigonometric and rational. Moreover most of the known models 
are linked to the last two cases only, i.e to the trigonometric r-matrix 

rt{X - A«) = ■ ■ ,] r f ^ C0S(A - //)(73 (8) (73 + (7+ (g) (7_ + (T_ (8) (7+) (4) 

I sm(A — /ij Vz / 
or to its g — > 1, A, — > limit given by the rational solution 

r,(A -ii)= .^^ ]_ ^^ P, where P ^ ^{I + a ■ a), (5) 

P being the permutation operator. The above remarkable observation has motivated us to 
conjecture that all integrable models satisfying the CYBE (2) must be derivable from an ances- 
tor model with their Lax operators obtained as various reductions of this single ancestor Lax 
operator and this should make the r-matrix, inherited from their ancestor, to be naturally the 
same for all these descendant models. In the next section we present such an ancestor model in 
the explicit form associated with the trigonometric r-matrix (4) , from which we will be able to 
generate a rich collection of integrable discrete and continuum models including inhomogeneous 
as well as hybrid systems, all satisfying the CYBE and sharing the same r-matrix (4) (or its 
rational hmit (5)). Note that from the CYBE (2) one can go to its global description 

{T^(A)®, T^ifi)} = [r(A - /x), T^(A) ® T^(/i)] (6) 

for the monodromy matrix 

T„(A) = MA)---MA)== (7) 

It is important to notice that (6) exhibits exactly the same form as its local relation (2), 
which reflects a deep underlying Hopf algebra structure, an important characteristic of all such 
integrable systems [15]. Defining now the transfer matrix as r(A) = irTjv(A) = aAr(A) -|- dAr(A) 
and taking the trace of (6) one can easily derive (since the rhs being the trace of a commutator 
is zero) the integrability condition (1) for the system. Therefore going backwards in the logical 
chain we can conclude that the nonlinear systems with its representative Lax operator and the 
r-matrix satisfying the CYBE (2) must be an integrable system. We shall see below that the 
relation (6) also carries important information for deriving action-angle variables and reflects 
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an universal property for all integrable systems sharing the same r-matrix and hence belonging 
to the same class. 

Note that in this algebraic approach we are not concerned about the usual Lax pair L, M 
and do not obtain the dynamical equation from the flatness condition involving them. We 
on the other hand take the Lax operator Lk{X) satisfying the CYBE as the representative of 
the integrable model and using it construct the monodromy matrix: T{X) — YikLkiX) and 
then the transfer matrix from its trace t(A) = trT{X). Expanding further the transfer matrix 
r(A) in spectral parameter A as described above, we derive the conserved quantities including 
the Hamiltonian H in the explicit form. The dynamical equation can now be obtained as the 
Hamilton equation ipt = {ip, H}, using the fundamental PB relations. 

At the lattice constant A ^ one may recover in some cases the corresponding field model: 
Lk{X) I + AC{x, A) + O(A^) with C{x, A) as the field Lax operator. Though the associated 
r-matrix remains the same, the CYBE gets deformed and the corresponding monodromy matrix 
T(A) at the infinite interval limit / = A^A — > oo satisfies also a bit different global CYBE [5]. 
For continuum models one can extract the conserved quantities more conveniently from the 
Lax operator using the Ricatti equation derived from the linear spectral problem. 



III. ANCESTOR MODELS ASSOCIATED WITH TRIGONOMETRIC AND RATIONAL r-MATRIX 

As mentioned, our generating scheme for integrable models is based on various reductions of a 
discrete ancestor Lax operator, which we propose to take in the following form [7] 

and demand it to satisfy the CYBE (2) with the trigonometric r-matrix (4). Sk appearing 
in (8) are the basic dynamical fields PB algebra of which as specified below is dictated by 
its integr ability and c^,a — 1,2 are a set of arbitrary parameters. The structure of the Lax 
operator (8) becomes clearer if we notice its possible decomposition, after an allowed gauge 
transformation by /i = e'^"^, L*(""^)(0 ^ /iL*(""^)(0/i~^ = ^L+ + ^-^L_, where L± are spectral 
parameter ^-free upper/lower triangular matrices. Note that the r-matrix (4) allows also a 
similar decomposition (after a similar gauge transformation): rt{^) — > ^r+ -|- (^)~-'^r_, ^ = 
e*^,77 = e^^ with r± being spectral- free upper/lower triangular matrices, which together with 
L± satisfy the FRT-type [16] PB algebra derivable from the CYBE [17]. The demand of 
integrability on (8) through the CYBE can be shown to be equivalent to the underlying general 
algebra 



{Si S^} = ±iSkiS^, Sf} = t^f{2aSl), with f{x) = (m+ sin(x) + M" cos(a;: 



(9) 
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where M"^ — ±|-\/±T(ci"c^ ± C1C2) are arbitrary parameters acting as central elements with 
trivial brackets with all others: {M='=, •} = and in general may also be site and time dependent. 
It is important to note that the underlying PB structure (9) is linked with a generalization of 
the well known quantum group algebra. For generating integrable systems from this ancestor 
model, we find first a realization of (9) in canonical variables {uk,pi} — S^i, in the form 

Sl = Uk, S+ = e-'P-g{u,), = g{uk)e^^\ (10) 

where ^ 

g{uk) = [n + sva.a{s-Uk){f{a{uk + s + l))}]^ , (11) 

sm Gt 

containing free parameters k, s and function f{x) as defined in (9). It should be remarked 
here that realization (10) usually assumes the complex conjugacy 5^ = (5'^)*, which however 
is not imposed by the integrabihty condition (9). Note that we have now lots of freedom for 
generating descendant models from the ancestor Lax operator (8) by using various reductions 
of (11) 

under different choices of the arbitrary parameters c's as well as k and s or its further 
realization in bosonic variables: {ipkjipf} — iSki in (10). Moreover we can multiply these Lax 
operators from left or right by (7a, a = 1,2,3 , since such transformations are allowed by the 
CYBE due to a symmetry of (4) and (5) as [r, aa <8) (Ja] = 0. 

We will demonstrate in the next section that a class of discrete integrable systems with 
nontrivial deformation parameter q, which may be interpreted as the relativistic parameter 
can be generated in a systematic way from the ancestor Lax operator (8). The nonrelativistic 
models on the other hand may be constructed in a similar way from the — > 1 limit of (8) 
given as 

-rat.{anc) f ^-^ _ ( ^^(A + s|) + cj, Sf. 

with c^'^,a = 0, 1 being arbitrary parameters. Here due to the corresponding limits of — > 
s, {cf} {c°'^}, M+ -m+, M- -am-, the PB algebra (9) reduces to 

{4> Sfe } = iSki{^m+sl + m"), {si, } = ±i6kis^, (13) 

where m'^ — c5c2, m" — c\c2+Cicl with {m^, •} = 0. Note that a Casimir operator commuting 
with all other generators of (13) may be constructed as 

S' = sl{m^sl + m-) + sts^ (14) 

and a realization of it (13) given by the generalized Holstein- Primakov traxisiormeition (HPT) 

sl^s-Nk, 4 = ^o(A^fe)V'fe, Sk^rk9o{Nk), go{Nk)^{m- + m+{2s-Nk))K Nk = rk%IJk 

(15) 
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in bosonic variables ipk, which in fact is the a ^ hmit of (10) and (11). We stress again 
that since the conjugacy of is not necessarily imposed by the integrability, il>,il>* in (15) 
in general may not be complex conjugates. Note that the ancestor model (12) represents the 
undeformed rational class and satisfies the CYBE with the rational r-matrix (5). (8) and (12) 
serving as the ancestor Lax operators for the discrete integrable models may also yield for some 
systems the corresponding field models with the Lax operator jC{x, A) . The associated r-matrix 
however would remain the same at the continuum limit, since it is a global nondynamical object 
independent of site indices. We shall see in sec. V that parameters c's in general can be space- 
time dependent and hence could induce inhomogeneity in the model preserving the constancy 
of the spectral parameter. 



IV. UNIFIED GENERATION OF DISCRETE INTEGRABLE MODELS 

Prom the ancestor models proposed we generate here integrable discrete models belonging to 
both trigonometric and rational class 



A. Relativistic models belonging to trigonometric class 

For constructing this class of models we start from the ancestor Lax operator (8) and look into 
its different realizations by choosing first the arbitrary parameters c's as constants. 

1. ) Discrete sine-Gordon model: Parameter choice cf = — = mA, with m as the constant 
mass. This gives M~ = 0, = — (mA)^, and reduces realization (10) correspondingly to yield 
from (8) (after multiplying it from right by — i(7i) the Lax operator 

/ Q{uk) e^^''^ , mAsin(A -I- awfc) \ o, , , .,9 , x . s 

Lk{\)^[ \ , 4. A / J' 9\uk)^l-{mAfcosa{2u,+ l). {16) 

y mAsm{X- auk), e '■^'''^ g[uk) j 

It is important to note that (16) yields exactly the Lax operator of the integrable discrete 
sine-Gordon model [18] and at the continuum limit A ^ 0, when 6^*^*=^ ^ 1 ± Aip^ and 
{uk,Pk) {u{x),p{x)), recovers clearly the field Lax operator 

Lfc(A) ~ l-|-A>C(a;, A), C{x,\) = ip{x)a3+m sm{\+au{x))a++m sm{\—au{x))a-, p{x) = u{x) 

(17) 

of the well known sine-Gordon model Uu — u^x + sin au — 0. Remarkably the PB algebra (9) 
in this case reduces to the classical limit of the celebrated quantum group [15] with its familiar 
relation = —i\2S'^\q. We will see in the next section that a more general choice for 

the parameters would lead to an inhomogeneous extension of this sine-Gordon model. 

2. ) Discrete Liouville model: Parameter choice cf — C2 — A, Ci — cf — 0. This gives 
— ±|-\/±T^^ and correspondingly reduces (10) to derive from the same ancestor (8) (after 
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multiplying it from right by ai) the Lax operator 



which represent the discrete Liouville model [19] and at its field limit (A — > 0) the Lax 
operator C — pa^ + e""(^(7+ + |cr-) of the well known Liouville equation: Uu — u^x — e"". 
Note that in this case (9) gives a novel PB algebra with exponentially deformed relation like 

L ' J 2 sin a 

It is intriguing to observe here that though the underlying PB structure and hence its 
realization giving the model are fixed by the choice of M^, the Lax operator (8) which depends 
directly on the parameters c's may take different forms for the same model. For example, in 
the present case with additional choice cj" 7^ would record the same values for M"^, but 
a different Liouville Lax operator [20]. This opens up therefore a promizing possibility for 
systematically obtaining different useful Lax operators for the same integrable model. 

3.) Relativistic Toda chain: Different sets of constant choices i) c+ = 1 , a = 1, 2, or ii) c~ = 
1 , a = 1, 2, or in) cf = ±1, or iv) cf = 1, with the rest of c's being zero, lead to M"^ = 0, 
reducing therefore (9) to the simple PB algebra 

{St, Sf} = 0, {Si S^] = ±i5,iSi, (19) 

and the realization (10) (after a canonical interchange of variables: u — > —ip,p — > —iu,) to the 
form 

S! = -tPk, Si = ae^^^ . (20) 

This generates interestingly from the same ancestor Lax operator (8) different forms of the 
discrete-time or relativistic Toda chain (RTC). For example, case iii) yields 

recovering the Lax operator found in [21], while iv) generates a different Lax operator [22] for 
the same model. More famous RTC model of Suris [23] however is obtained in this approach 
after performing a twisting transformation with twisting parameter taken as ±a ( of these 
equivalent cases we consider here only —a, for definiteness) , which deforms the r^- matrix (4) 
by adding a constant matrix fl to it [21]: 

ft^fQ — rt — where Q = i((T3 (8) 7 — / (g) (73). (22) 

As a result the form of the ancestor Lax operator (8) also gets changed with its elements 
transforming as 

- c^e-^S Si Si = e^^«^^ e^'^^'k (23) 
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Implementing the corresponding changes in (21) and using the same reahzation (20) for the 
variables (23) we obtain now the explicit form of the Lax operator 

generating that of the Suris RTC [23] . 

4. ) Discrete derivative NLS: Parameter choice as constants: 

iA 

c+ = = 1, c7 = -igA, C2 = — , giving M+ = 2AsinQ;, M = 2iAcosQ; (25) 

gives (10) a q-bosonic realization as S'/^ — —Qk, — Ql, — —Nk, with a PB algebra 
induced from (9) as 

{Qk, Ni} = zSkiQ, {Qk, Qn = iSki^^^^^^^^^, (26) 

cos a 

which clearly reduces to the standard bosons ■0^, V'fe o; — > 0. It is worth noting that this new 
q-bosonic model as realized from the ancestor Lax operator (8) (after introducing A) would 
give 

which represents an exact lattice version of the DNLS equation [24] . When expressed through 
bosonic field: Q = N = A|V'P, (27) yields at the continuum limit A ^ the 

field Lax operator 

^ W = -ile - m')<J3 + e(V'V+ + V^a-) (28) 

of the well known Chen-Lee-Liu DNLS equation [25]: iipt = ipxx — 4:i\ip\'^ipx- 

5. ) Ablowitz-Ladik model: This model involving also another form of g-boson is possible to 
generate in our scheme, though it needs twisting transformation as in the Suris RTC mentioned 
above and is associated with the same twisted rj^-matrix (22) and the twisted ancestor Lax 
operator with the change (23). Now the the parameter choice Ci = C2 = with Ci = C2 = 1 
giving M"^ = |V=tl (compare with the Liouville case) together with the twisting removes 
dynamical variables from the diagonal elements of the twisted Lax operator as well as modifies 
the Poisson algebra of the transformed variables as derivable from (9). Therefore naming 
bk — 2sinQ;5'^ we get this modified PB relation as {bk, = iSki{l — blbk) , confirming the 
basic variables of the Ablowitz-Ladik model as a type of g-boson with its Lax operator as 



L,{0 =15' \. (29) 
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related to (22). We will see later how space-time dependent parameters c's give inhomogeneous 
extensions of this model. Note that another intriguing possibility of generalizing this model 
arises if we simply consider Ci" 7^ in the above construction. It is not dificult to see that, 
realizing — — ln(l — blbk) this would generate an extra term ^Ci{l — blbk)~'^^°' in the upper 
diagonal element of the Ablowitz-Ladik Lax operator (29). Its consequence in the dynamical 
equation would be an interesting problem to study. 



B. Nonrelativistic models belonging to rational class 

Deformation parameter q = e*°, as we have seen in the above models, serves as the relativistic 
or the deformed bosonic parameter. We consider now the undeformed limit g — 1 or a — , 
when as explained already, the r-matrix reduces to its rational form (5) and the ancestor Lax 
operator is converted to (12) with the underlying PB algebra (13). 

We find that the integrable models belonging to this rational class are mostly nonrelativistic 
models, which can be generated in a similar way from the rational ancestor model (12) with 
different constant choices for parameters c°'^, a = 1, 2 involved in it. 

6. ) Landau- Lifshitz equation (LLE) Parameter choice — l,c\ — —c\ — —I compatible 
with m+ = l,m~ = 0, reduces (13) to the classical 5/2 spin algebra {s^,sf} = i5kie°'^'^ with 
spin: = (s|)^ + s^s'^ = s| as the Casimir operator reduced from (14). The ancestor Lax 
operator (12) simplifies ( ignoring an irrelevant multiplicative factor) to 

LkiX) ^ I + j^^Sk ■ B (30) 

representing a discrete version of the LLE. At the continuum limit A ^ putting Sk As{x) 
one gets from the Casimir: ^{x) = 1 and from the Lax operator Lk{X) — >■ / + A£{X), JC{X) — 
^s(x) ■ a. that for the well known LLE [26]. 

7. ) Discrete NLS model : For the same s/2 spin algebra transformation (15) yields the 
standard HPT with go{\^k\) = (2s — A|-?/'fcp)5 (considering ipkii^t ^'^ t)C complex conjugates 
and scaling them by ^J~K). This realization by considering parameter / = leads (12) to the 
Lax operator of exactly integrable discrete NLS model [18] given by 



Lk{X)=\ ''r.: .T' "^".T"; I- (31) 



At the field limit: A ^ 0, (31) yields (after multiplying it from left by a^A and considering 
s = 1/ A) the familiar form of the Lax operator 

£(A) = Xa^ + V2{ij*a+ - V'cr") (32) 

for the NLS field equation iipt = ipxx + 2\ip\'^ip. 
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8.) Simple lattice NLS: On the other hand a complementary choice m'^ — 0,m~ = 1, 
giving 5'o(-A^fc) — 1 converts (15) directly to the realization = ■0^, — i^l, ^l — s — iplipk in 
bosonic field : {■0fe, V'f} — '^^ki- Now a compatible choice of parameters: — cl — 1, C2 — c\ — 
together with this bosonic realization generates from the ancestor (12) the Lax operator 



which may be associated with another simple lattice NLS model proposed in [27] and as 
noted there ip,ip* may not be complex conjugates at the discrete level. At the continuum limit 
we recover again the same field Lax operator (32) for the NLS equation and regain also the 
complex-conjugacy of the fields (see for details [27]). 

9.) Nonrelativistic Toda chain: Note that the trivial choice = yields from (13) again 
the same algebra (19) and therefore we may take the same realization of it as found before. 
However the rational form of ancestor model (12) generates now simpler Lax operator 



of the nonrelativistic Toda chain associated with the rational r-matrix and described by the 
Hamiltonian: -H" = Efc IpI + e("'=-'"=+i) . 

Thus we have demonstrated that discrete and continuum intcgrable models can be obtained 
in a unified way from the ancestor Lax operator (8) (or its rational limit (12)) by choosing 
different sets of constant values for the parameters c's involved in the ancestor model and by 
using different realizations of the underlying PB algebra. In the next section we find how a 
more general choice of c's can generate further the inhomogeneous extensions of these integrable 
models. 

We find also a convincing answer to an important question raised above asking why differ- 
ent integrable systems with varied Lax operator solutions should have the same r-matrix, by 
discovering that all these models are basically obtainable from the same ancestor model (8) 
associated with the trigonometric r-matrix (4). These descendant models, whose explicit Lax 
operators we derive here satisfy the GYBE (2) inheriting and sharing the same r-matrix (4). 
We will see in sec. VI that this significant fact induces a universality among these seemingly 
diverse systems by defining their action-angle variables in the same way. 

V. INTEGRABLE INHOMOGENEOUS AND HYBRID MODELS WITH ISOSPECTRAL FLOW 

We have seen how by fixing the values of certain parameters we could generate a wide spectrum 
of integrable models belonging to the trigonometric and rational class. We focus here on some 
promising possibilities to generalize this procedure for constructing novel integrable families of 
inhomogeneous and hybrid models. 




(33) 




(34) 
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A. Inhomogeneous models 

Returning again to the ancestor model (8) we may notice that the parameters c^,a = 1,2 en- 
tering in it (similarly, c|^'^, a = 1, 2 in its rational limit (12)) act like external parameters having 
trivial PB with all basic variables in their local algebra and therefore, apart from constants as 
earlier they may be considered in general as site (time) dependent arbitrary functions. As a re- 
sult in (10), (11) in turn also become functions M^{t), k = 1,2, . . . , N (and similarly m^{t) 
in (15)) and lead to new integrable descendant models, which are inhomogeneous extensions of 
the discrete and continuum models constructed above. However this integrable family of inho- 
mogeneous models is obtained in our scheme by keeping the usual isospectral flow. Moreover, 
such constancy of spectral parameters (except some trivial transformations like shifting etc.) is 
essential in this algebraic formalism for satisfying the CYBE (6) with spectral dependent global 
and nondynamical r(A — /x)-matrix. It is important also to notice that the inhomogeneity is 
introduced here through a set of different independent parameters: c^j. (or c^^ ) with a = 1, 2 
and therefore it may not be always possible to absorb them in the single spectral parameter, 
even by declaring it to be nonisospectral. Therefore we see that, contrary to the standard 
approach the inhomogeneous models can not be described in general as nonisospectral flow, at 
least those that belong to the present family. Moreover isospectrahty is a necessary criterion 
for the CYBE solution, as explained already. 

10. ) Variable mass sine-Gordon model: 

The construction is parallel to that of the constant mass sine-Gordon model obtained above, 
where in place of constants we choose now the parameters as four different variable mass: 
cf = mfi^(t)A, cf = mfi^{t)A. This would generate from the ancestor Lax operator (8) and 
realization (10) a general form of a new inhomogeneous sine-Gordon model, which is integrable 
and satisfies the GYBE associated to the trigonometric r-matrix (4). Particular choices of the 
inhomogeneities would yield naturally different forms of the variable mass sine-Gordon model, 
discrete as well as continuum, which seem to have been never considered before. 

For a demonstration we take up the simplest case when all mass parameters coincide: 
m^^.(t) = mk{t). This variable mass discrete sine-Gordon model can be described again by 
the same form (16) by replacing constant m by a variable mk{t). At the continuum limit this 
would correspond to a sine-Gordon field model with variable mass m{x, t). If the mass parame- 
ter is assumed to be factorized: m{x,t) = mQ{t)mi{x), by introducing a new coordinate system 
through nonlinear transformation {t,x) — > (T,X), T = /*mo(t')rft', X = f^mi(x)dx', the 
Hamiltonian of the model can be written formally again as the standard sine-Gordon model 
with unit mass. Nevertheless we notice that even in a further simplified case with mo{t) — 1, 
the sohton solutions in the original system might have quite interesting character depending 
on the form of the variable mass mi{x) (see Fig. 1). 

11. ) Inhomogeneous NLS model: Since this model belongs to the rational class, in accordance 
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with our strategy we start with the ancestor Lax operator (12) and consider the parameters 
involved in it and in reahzation (15) to be site and time dependent functions: c°-^^(i), a — 1,2. 
With all of them different we naturally get the general inhomogeneous discrete NLS model, 
which retains its integrability and contrary to the standard approach also its isospectrality, as 
explained already. For constructing the corresponding field model we take the parameters in 
the form c?^ = gu, c^k = -fl'2fe, = ^ + fik, = + f2k, s = A and at the limit A ^ 
obtain the Lax operator 

C{\)= I ^^^"^'^^ ^\ ), where A„(x, i) = A^„(x, i) + /„(x, i), a = 1,2 (35) 
\^ -Q -A2{x,t) J 

and Q — ipgoix, t) with go{x, t) — {gi{x, t)+g2{x, t))^, representing an inhomogeneous NLS field 
model with inhomogeneities introduced by the independent functions ga{x, t), fa{x, t),a — 1,2. 
It may be stressed again that here the spectral parameter A is strictly constant and when Ai ^ 
A2, all inhomogeneous parameters apparently can not be absorbed in this single parameter. 
It is challenging to derive the explicit form of this integrable variable coefficient general NLS 
equation, associated with the rational r-matrix (5). For showing that the Lax operator of 
many known inhomogeneous NLS equations can actually be derived from (35), we consider the 
particular situation gi — g2 = g{x, t), fi — f2 = f{x, t) and rewrite (35) as 

C{X) = A{x,t)as + Q*a+-Qa_, A{x,t) = Xg{x,t) + f{x,t). (36) 

It is remarkable that from this single operator we recover at g — 1, f — at the Lax operator 
of [9] , at = |, / = ^ that of [10] and similarly at ^ = T{t), f = ^xT{t) that of [11]. Note 
however that the actual form of the equations depend also on the time evolution operator M, 
which is likely to be different in our approach from the known ones, since in our construction the 
fundamental canonical PB structure is always preserved. Therefore it would be a challenging 
problem to derive these new integrable inhomogeneous NLS equations explicitly from their 
Hamiltonian using the canonical PB. 

12. ) Gaudin model: It is intriguing that by just by considering the parameter / in the Lax 
operator (30) for the discrete LLE to be site dependent: I ^ Ik, one recovers the Lax operator 
Lk{X) = I + j^^k ■ <? for the celebrated Gaudin model, given by the integrable Hamiltonians 
Hk = Yjh^k j^Tii^k ■ si), k — 1,2, . . . ,N [28]. This model is associated also with the rational 
matrix (5). 

13. ) Inhomogeneous relativistic and nonrelativistic Toda chains: It is not difficult to see 
that by repeating the construction of the Toda chains but taking the parameters to be noncon- 
stants we get integrable inhomogeneous extensions of such models. For example, considering 
in ancestor Lax operator (12) the parameters to be cf = /^(t), cf = 0, but using the same 
realization as for the original relativistic Toda chain, we get an extension of its Lax operator 
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(21) to include inhomogeneity through arbitrary functions fk'it)' 




which therefore would represent a new integrable family of inhomogeneous relativistic Toda 
chain. We will not present here its explicit form. 

At a — s> this family of relativistic models would go to its nonrelativistic limit represented 
by the Lax operator 



which is an obvious extension of (34) (by introducing the inhomogeneous parameter g2k = 

Ik 

and normalizing it by c^^)- Without defining any time evolution operator M, we can directly 
construct from (37) the explicit form of the Hamiltonian through the conserved quantity as 
H — Cn-1 and derive the Hamilton equations using the canonical PB between Uk,pi, yielding 
Uk — Pk + g2k and hence the inhomogeneous Toda chain equation as 



with arbitrary parameters gi{k) = (cij^c^j^^i)^^ and g2{k). Different choices of these parameters 
would generate from (38) different inhomogeneous Toda chains. For example the particular 
choices: gi{k) = k,g2{k) = a^t and gi{k) = 4A;^ + l,g2{k) = kt derives the Toda chains found 
in [12], though in contrast we recover this result in a completely isospectral way. 

14.) Inhomogeneous Ablowitz-Ladik model: It is easy to notice again that if instead of 
constants as in the original model, we choose the parameters through arbitrary function r(t) 
as Ci = {c2)~^ — e^^^\ keeping the same trivial choice for cf = C2 = 0, we generate from (8) 
the Lax operator 



Remarkably, in spite of our isospectral approach, (39) recovers exactly the Lax operator of [13] 
for arbitrary T{t) and that of [14] for T{t) = at, representing known inhomogeneous Ablowitz- 
Ladik models. 

In a similar way by generalizing the constant parameters to inhomogeneous functions one 
can generate systematically inhomogeneous extensions of other integrable models constructed 
here. Note again that all such extensions retain the integrability of the system as well as the 
isospectrality and the same r-matrix solution. 




(37) 



df 



Uk = ^i(A;)e"('=-^)-'^('=) - gi{k + i)e«('=)-'^('=+i) + g2{t) + boundary terms 



(38) 




(39) 
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B. Integrable hybrid models 

Our scheme for generating different integrable models from an ancestor model sharing the same 
r-matrix opens up a possibility of constructing new families of integrable models by hybridizing 
these descended models. 

Such constructions can be of two types. The first type of hybrid models may be constructed 
by using different descendant Lax operators obtained directly from (8) (or alternatively from 
(12)) as its different but consistent reductions and realizations at different lattice sites. Since all 
representative Lax operators of these constituent models: L'l^''\X), with d{k) denoting different 
members of the same descendant class inserted at sites k, should share the same r-matrix, the 
monodromy matrix of this hybrid model: T'^'^^(A) = Yld(^k),k Lf\X) must satisfy the global 
GYBE (6) and represent therefore an integrable system with the set of conserved quantities 
including the Hamiltonian obtainable as usual through expansion of t'^^^(A) = tr{T^^^{X)) in 
the spectral parameter. One can generate in this way some exotic hybrid models by combining 
for example, sine-Gordon and Liouville models, different types of relativistic Toda chain or 
discrete NLS model etc. constructed above. These hybrid models presumably would show 
different dynamics at different domains in the coordinate space. It is encouraging to note that 
very recently such models have received well deserved attention, though only at the continuum 
level [8]. We hope that the present idea, based on discrete approach and r-matrix formalism 
would prove to be promizing and fruitful for analyzing such hybrid integrable models. 

A second type of hybrid models may be constructed by considering different representation 
of the Lax operator for different components of the field and inserting their direct product 
at the same lattice site. As a result one can build new mult i- component generalization of 
a scalar model through the fused Lax operator: L^"^^ = YlrnL^k^\ where each entries in the 
product would represent individual components. Note that unlike the vector generalization, 
which needs also enlarged matrix realization for the Lax operator, our multi-component hybrid 
models would yield only 2x2 matrix Lax operators. For elaborating this idea we present the 
detail construction of an integrable hierarchy of two-component DNLS model. 

15.) Integrable hierarchy of two- component DNLS : Note that in constructing the discrete 
DNLS model in sec. IV, the values of M"^ fixed by (25) actually determined the underlying 
algebra as well as the required realization. It is however crucial to notice now that interchang- 
ing the parameters cf ^ 62 would not change the values of M"^ and therefore would lead to 
the same algebra and its realization, but result to a complementary form for the Lax oper- 
ator, though representing the same DNLS model. In our construction of the two-component 
model with fields = 1,2 having PB relations {'4'k'\i'*i''^'^} = "i^isi^kh we take c^, a = 1, 2 

as in (25) for building the Lax operator l[,^^(^*^^)) as (27) for the first component. However, 
for the corresponding construction of related to the second component we take the 

complementary choice by considering cf . The fused Lax operator taking the form 
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= L^^^('0(-'^))L^^^('0(^^) represents now a new discrete multi-component DNLS satis- 
fying the CYBE (2) with the same r- matrix (4). At the continuum hmit A — > 0, repeating the 
construction for the scalar DNLS model (28), it is easy to see that the field Lax operator of this 
two-component model is given simply as a linear superposition -0^^^) = 

with the explicit form 

^^^'\^^'^) = i\iY^-e) + wi-MW3+ m+-^r2W^ + {^^i+-^^2)a-, e = e^' (40) 

It is interesting to show by direct construction that this Lax operator generates an integrable 
hierarchy of multi-component DNLS model through the expansion Inr(A) = Z^^o^in-^^^" 
with H2 — C2 + C_2 introducing a new two-component generalization of the Chen-Lee-Liu 
DNLS equation. We are not giving here its explicit form, which can be worked out with little 
patience. The higher conserved quantities will yield higher order equations. Some similar class 
of discrete matrix and multi-component DNLS models was proposed recently [30] . 

16. ) Massive Thirring model: It is remarkable that (40) constructed through a novel hy- 
bridization from our ancestor model coincides also with the Lax operator of the bosonic massive 
Thirring model [29]. Hamiltonian of this relativistic model may be given through the same con- 
served quantities constructed for the above model in the form Hi — Ci + C_i. 

17. ) Integrable discrete self trapping model: A discrete self trapping model with two-bosonic 
modes ■^^"^ ■0*^"), a = 1, 2 given by the Hamiltonian 

2 



1 ^(5^ _ 7V(»))2 + (^*(l)^(2) + ^*(2)^(1)^) 

2 ^ 



was studied in [31]. This integrable model is associated with a Lax operator, which may be 
constructed by fusing two operators as L(X) = L^-^\X) L^'^^X) , where L^"'\X) are given by the 
same Lax operator (33) for each of the modes a = 1, 2. 

An interesting line of investigation would be to apply this hybridization method for con- 
structing possible multi-component extensions of other models like relativistic Toda chain, 
Abowitz-Ladik, Liouville model, LLE, Gaudin model etc. The linear superposition of Lax op- 
erators for building new nonlinear integrable systems, as revealed here, seems to be a promizing 
idea worth pursuing. 



VI. UNIVERSAL PROPERTIES OF INTEGRABLE DESCENDANT MODELS 

We have seen that diverse forms of integrable models: discrete and continuum, homogeneous 
and inhomogeneous, multi-component and hybrid models can be generated in a systematic 
way in our ancestor model scheme. Among this diversity however we find also an unexpected 
universality. Indeed, as we have found, a wide range of models, namely sine-Gordon, Liouville, 
DNLS, relativistic Toda chain etc. including their discrete, inhomogeneous and hybrid variants 
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belong to the trigonometric class, while models like NLS, Toda chain etc. and their related 
discrete and inhomogeneous extensions are in the same rational class, which being in fact the 
undeformed — > 1 or the nonrelativistic limit of the former class. 

The crucial observation is that the diversity of all descendant models belonging to the 
same class seems to disappear at the global level allowing their description through a universal 
action-angle variable. The reason for this is very simple. Though these models differ widely 
at their local level having different forms of the Lax operator, their monodromy matrix 7jv(A) 
(7) satisfies the same global relation (6) with the same r-matrix, which is inherited from their 
ancestor model and shared by all of them. 

As a result, for all models belonging for example to the trigonometric class, the PB relations 
should be given by the same structure constants expressed through the elements of the rt-matrix 
(4). For the twisted models, e.g. Suris RTC and Ablowitz-Ladik model the structure constants 
should similarly be given by the twisted r^-matrix (22). In the same way all models from 
the rational class should have analogous property expressed through the elements of rational 
rj.-matrix (5). However, while the action variables are constants in time, the time evolution 
of angle variables depends on the definition of the Hamiltonians through conserved quantities, 
which usually differs for different models. 

Such differences also bear some additional imprint at the continuum limit, when the mon- 
odromy matrix is defined as T(A) = limN^oo {X)Tfyf{X)L^ {X) and the corresponding 
CYBE is modified as [5] {T(A)®, T(/i)} = r+(A - /x)T(A) ® T(/i) - T(A) ® T(/i)r_(A - //), 
where r±(A - /x) = limTv^oo L^^W ® L^^{^^)r{X - /i)L±^(A) ® L^^i^i). Therefore though the 
action-angle description for such models of the same class are again basically the same, the 
influence of the individual models also enters now due to the appearance of r± modified by the 
asymptotic forms L^o^X) of the individual Lax operators. Nevertheless it is startling to check 
that the canonical action-angle variables for widely different field models hke DNLS [24] and 
the sine-Gordon [1] are defined exactly in the same way: p(^) = 2^^^ l'^(OI) ^(0 — ^^sKO 
for the continuum modes with PB {q{C),p{r])} — — r]), ^ > and = ^ In^fe, Qk — Inftfe 
for the discrete set with {qk,Pi} — Ski ^-nd similarly for their conjugates qk,Pi- 

Therefore we may conclude that all integrable models presented here may be described 
universally through the ancestor Lax operator (8) and the r^-matrix (4) (or twisted rn) or 
similarly by the 5 — > 1 limit as the ancestor model (12) and the r^-matrix (5), where the 
global relations like action-angle variables are determined by the r-matrix elements alone. The 
individuality of the models may be refiected only in the definition of their Hamiltonians through 
conserved quantities and for continuum models, additionally in the limiting forms of their Lax 
operators. 
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VII. CONCLUDING REMARKS 

Wc have presented here a unifying scheme based on PB algebra for systematically generating 
a large class of integrable discrete and continuum models from a single ancestor model. Such 
models include well known and new integrable systems as well as inhomogeneous models. Based 
on our construction we conclude that more general and logical approach for inhomogeneous 
integrable models , at least for models with nondynamical r-matrix, would be to describe them 
as isospectral flow in inhomogeneous external flelds. As another fruitful application of the 
present scheme we have proposed a simple method for constructing new families of integrable 
hybrid models by fusing different types of descendant models. In spite of the vastly diverse form 
of these models their common ancestor and common r-matrix reveal an inherent universality 
in their description through action-angle variables. 

We strongly hope that the algebraic approach linked with the quantum group structure 
formulated here for generating classical integrable discrete as well as field models, though a 
bit uncommon in the community working in classical integrability, would prove to be much 
powerful due to its systematic and algorithmic nature. Similarly the novel ideas of construction 
introduced here, like generating integrable hybrid and multi-component models and creating 
integrable inhomogeneity in isospectral fiow, are expected to be equally promizing. 

References 

[1] S. Novikov, V. Manakov, L. Pitaevskii and V. Zakharov, Theory of Solitons (Plenum, N.Y., 
1984) 

[2] L. J. Mason and G. A. J. Sparling, Phys. Lett. A 137 29 (1989) ; S. Chakravarty, M. J. 
Ablowitz and P. A. Clarkson, Phys. Rev. Lett. 65 1085 (1990) 

[3] V. I. Zakharov and L. A. Takhtajan, Teor. Mat. Fiz. 38 (1976) 26 
Anjan Kundu, J. Math. Phys. 25 (1984) 3433 

[4] Symmetries and Integrability of Difference Equations, (ed. P. A. Clarkson and F. W. Ni- 
jhoff, Cambridge Univ, Press, 2001) 

[5] L. D. Faddeev, in Recent Advances in Field Theory and Statistical Mechanics (ed. J. Zuber 
and R. Stora, North-Holland, Amsterdam, 1984), p. 561 

[6] Anjan Kundu and B. Basumallick Mod. Phys. Lett. A 7, 61 (1992) 

[7] Anjan Kundu, Phys. Rev. Lett. 82, 3936 (1999) 



18 



[8] P. Bowcock, E. Corrigan and C. Zambon, Classically integrable field theories with defects, 
ArXiv:hep-th/0305022, May 6, 2003 



[9 

[lo; 
[11 

[12 

[is; 

[14 
[15 

[16; 

[17; 
[is; 

[19 

[20 
[21 

[22 
[23; 
[24; 
[25 
[26 
[27; 



H. H. Chen and C. S. Liu, Phys. Rev. Lett. 37, 693 (1976) 

R. Radha and M. Lakshmanan, J. Phys. A 28, 6977 (1995) 

R. Balakrishnan, Phys. Rev. A 32, 1144 (1985) 

D. Levi and O. Ragnisco, J. Phys. A 24, 1729 (1991) 

V. Konotop, O. Chubykalo and L. Vaazquez, Phys. Rev. E 48, 563 (1993) 

R. Scharf and A. R. Bishop, Phys. Rev. A 43, 6535 (1991) 

V. G. Drinfeld, Quantum Groups in Proc. Int. Cong. Math. Berkeley, 1986 (Americal Math. 
Soc, Providence, 1986) 1 , 798. 

N. Yu. Reshetikhin, L. A. Takhtajan and L. D. Faddeev, Algebra and Analysis 1, 178 
(1989) 

Anjan Kundu, Theor. Math. Phys. 99, 428 (1994) 

A G Izergin and V E Korepin, Nucl. Phys. B 205 [FS 5], 401 (1982) 

L. D. Faddeev and L. A. Takhtajan Lect. Notes Phys. 246 (ed. H. de Vega et al, Sringer, 
Berlin, 1986), p. 166 

L.D. Faddeev, O. Tirkkonen, Nucl. Phys. B 453, 647 (1995) 
Anjan Kundu, Phys. Lett. A 190, 73 (1994) 
R. Inoue and K. Hikami, J. Phys. Soc. Jpn, 67, 87 (1998) 
Yu. B. Suris, Phys. Lett. A 145, 113 (1990) 

Anjan Kundu and B. Basumallick, J. Math. Phys. 34, 1252 (1993) 
H. H. Chen, Y. C. Lee and C. S. Liu, Phys. Scr. 20, 490 (1979) 

H. C. Fogedby J. Phys. A 13, 1467 (1980) 

I. Merola, O. Ragnisco and Tu. G. Zhang, Inverse Prob. 10, 1315 (1994) 
Anjan Kundu and O. Ragnisco, J. Phys. A 27, 6335 (1994) 



19 



[28] E. K. Sklyanin, J. Sov. Math. 47, 2473 (1989) 

A. Hone, V. Kuznetsov and O. Ragnisco, Backlund transformation for the sl(2) Gaudin 
magnet, eprint: nlin.SI/0007041 (2000) 

[29] P. Kulish and E. K. Sklyanin, Lect. Notes in Phys. (ed. J. Hietarinta et al. Springer, Berlin, 
1982) 151, p. 61. 

[30] T. Tsuchida, J. Phys. A 35 7827 (2002) 

[31] V. Enol'skii, M. Salerno, N. Rostov and A. Scott, Phys. Scripta 43, 229 (1991) 



20 




Fig 1: How the kink solution (for mi{x) = 1) deforms in variable mass sine-Gordon model 
depending on the mass parameter mi (x) = x and mi (x) = , respectively. 
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